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We have developed a formalism to study non-adiabatic, non-radial oscillations of non-rotating 
compact stars in the frequency domain, including the effects of thermal diffusion in the framework of 
general relativistic perturbation theory. When a general equation of state depending on temperature 
is used, the perturbations of the fluid result in heat flux which is coupled with the spacetime geometry 
through the Einstein field equations. Our results show that the frequency of the first pressure (p) 
and gravity (g) oscillation modes is significantly affected by thermal diffusion, while that of the 
fundamental (/) mode is basically unaltered due to the global nature of that oscillation. The 
damping time of the oscillations is generally much smaller than in the adiabatic case (more than 
two orders of magnitude for the p— and g— modes) reflecting the effect of thermal dissipation. Both 
the isothermal and adiabatic limits are recovered in our treatment and we study in more detail the 
intermediate regime. Our formalism finds its natural astrophysical application in the study of the 
oscillation properties of newly born neutron stars, neutron stars with a deconfined quark core phase, 
or strange stars which are all promising sources of gravitational waves with frequencies in the band 
of the first generation and advanced ground-based interferometric detectors. 



I. INTRODUCTION 



The theory of stellar oscillations has been a fundamental tool in the study of stellar interiors and stellar properties 
during decades. For Newtonian stars the theory is well established and observationally tested. In some areas such 
as helioseismology, very high precision measurements of the normal oscillation frequencies allow for a detailed under- 
standing of the properties of the solar interior (see [1] for a review) . It is usually assumed that stellar pulsations are 
adiabatic because the thermal relaxation timescale in stellar interiors is orders of magnitude larger than the pulsation 
periods. However, in some particular situations, energy transfer in the external regions of stars is fast enough to 
affect the pulsation properties, and some work has been devoted to study, for example, non-adiabatic oscillations of 
white dwarfs [2] or the coupling between the different non-linear modes in non-adiabatic situations [3]. Concerning 
relativistic stars, the study of pulsation properties of neutron stars or compact objects such as strange stars or hybrid 
stars has become more popular in the last decade (see e.g. the reviews in Ref. [4-6]), mainly due to the expecta- 
tions of detecting the gravitational emission from pulsating nearby compact stars with the current or next generation 
ground-based interferometric detectors (LIGO, VIRGO, GEO600, TAMA). The theory of non-adiabatic relativistic 
stellar pulsations is not as well developed as its Newtonian cousin, due the higher complexity of the formalisms, 
but also due to the fact that most attention has been paid to the study of old neutron stars, in which the thermal 
conductivity is too small to have visible non-adiabatic effects. Nevertheless, there might be situations in which this is 
not entirely true. For example, in a newly born neutron star the thermal structure is determined by neutrino diffusion 
[7,8], instead of electron conduction as in old neutron stars [9], and the effects of non-adiabaticity are likely to be 
relevant in the outer layers. Another interesting possibility is the existence of deconfined quark matter in neutron 
star cores, or in the form of strange stars. There seems to be common agreement in that, if strange matter exists, 
it has to be in a color superconducting phase. It has been recently pointed out [10] that the thermal conductivity 
of such exotic matter is many orders of magnitude larger than in normal neutron star matter, and for sufficiently 
high temperatures the timescale for thermal relaxation can be as short as 10 -4 s., comparable with typical oscilla- 
tion periods of compact objects. If this or other exotic scenarios (hyperons, kaon condensates) happen to be true, 
our common belief that compact star pulsations are adiabatic must be modified, and some care must be taken to 
understand how non-adiabatic effects change the oscillation properties and therefore the predicted gravitational wave 
signals of pulsating compact stars. 

With this motivation, we have derived a formalism that includes the effects of heat transfer and chemical diffusion 
(important in proto-neutron stars, where lepton diffusion is the driving force of thermodynamical changes) in a 
relativistic analysis of stellar perturbations. Some work in this line has been done in the past for radial oscillations 



1 



[11]. We have considered the case of non-radial oscillations, since this is the case of interest for gravitational wave 
emission, by extending in a simple way the formalism of Lindblom and Detweiler [12,13], and complementing the 
system of equations in the frequency domain with the additional equations for thermal or chemical diffusion. In the 
present study we focused on the effects of heat transfer and chemical diffusion, neglecting the rotation of the star. 

The paper is structured as follows. In Sec. II we derive the equations of non-adiabatic stellar perturbations in 
general relativity. In Sec. Ill we describe the equation of state we have used and we dehne the thermodynamical 
quantities we need in our derivation. The additional equations (transport of energy) that close the system are discussed 
in section IV. Section V is devoted to the numerical implementation of the complete set of equations. In Sec. VI we 
discuss in detail the results of the numerical integrations and in Sec. VII we draw the main conclusions and comment 
on possible future extensions of this work. 

II. DERIVATION OF THE EQUATIONS 

The stress-energy tensor of a non-perfect fluid 1 , including heat flux but without viscosity, has the general form [14] 

T a (3 = {p + P)u a U[3 + pg afi + u a q ( j + u a qp (2.1) 

where p is the energy density, p is the pressure, u a is the matter four-velocity, and q a is the heat flux which satisfies 
u a q a — . In addition, we will also consider the conservation equation for the baryon number density n (equation of 
continuity) 

(nu a ). a = 0. (2.2) 



A. Background configuration 

Hereafter, we will neglect the heat flux in the background configuration. This is justified if the background is assumed 
to be in thermal equilibrium. Even in the case that thermal or chemical gradients are present, the assumption of 
stationary background is valid if the timescale for global thermodynamical changes is much larger than the timescale 
of variation of the perturbations. For example, in newly born neutron stars, structural changes happen on timescales 
of the order of 0.1-1 second, while we are interested in oscillations of periods of the order of milliseconds. Therefore, 
labeling background quantities by the superscript (°) , in the following we assume that 

= (2.3) 

such that, our background is a perfect fluid, spherically symmetric star, and is described by the well known TOV 
equations: 

g(0) 



u 



diag(- e ^\ e A M,r 2 7afc ) (2.4) 

(0) ^ = (e-' y / 2 ,0,0,0) (2.5) 

A' = -^e A (M-47r^°V 3 ) (2.6) 

v ' = ;^e A (M + 47r/V) (2.7) 

P (0) ' = -i(p (0) +P (0 V,, (2.8) 

where we denote with a prime d/dr, and we have defined 

lab = diag(l,sin 2 tf) (2.9) 

(here and in the following, greek indexes \i = 0, . . . , 3 run on spacetime, latin indexes i = 1, . . . , 3 run on the spatial 
subspace, and latin indexes a — 9,<j> run on the sphere). 



throughout the paper, we use units in which the Newton gravitational constant, the speed of light, and the Boltzmann 
constant are equal to one, G = c = &b = 1. 
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B. Perturbations 



The equations of stellar perturbations in the perfect fluid case have been derived by many authors in different 
formalisms [15,16]. In this paper, we follow the notation of Lindblom & Detweiler [12,13], (LD hereafter) and we will 
try to take the parallelism between their and our equations as far as possible. The perturbed stress-energy tensor 
has the form 



ST^ = (Sp + &p)ufu<® + (p (0) + p^)(5u^ + uf8u v ) 
+5pg ( °J +P (0) Sg MV + (Sq^ + u^6q v ) . 



(2.10) 



Following the conventions in [12,13], we expand in spherical harmonics the metric perturbations with polar symmetry 
(we will not discuss perturbations with axial symmetry, which are not related with fluid oscillations in non-rotating 
stars) as 



(2.11) 
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The Lagrangian displacements are expanded as 

r = - 

e 



-A/2 



-W lm Y lm r l e 



iuit 



and they are related to the four-velocity perturbation by 
so that its expansion in harmonics is 



(* = 1,2,3) 
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We also expand the perturbed heat flux in harmonics as follows 

5q a = Ke-»l 2 (0,Qu m (ry- 1 Y lm ^^),Q 2lm (r)r l Y^(d^)) , 
(k thermal conductivity) and the Lagrangian perturbations of the pressure and the energy density as 

Ap = Sp + pM? = Ap lm Y lm r l e iut 



(2.12) 



(2.13) 



(2.14) 



(2.15) 



Ap = Sp + p^e = Ap lm Y lm r l e 



I ml Atot 



(2.16) 



(the same holds for the Lagrangian perturbations of the other thermodynamical quantities) . 

Finally, following and generalizing LD, we will use a set of variables (Xi m , E[ m , £; m ) related to the Lagrangian 
perturbations of the pressure, energy density, and entropy, defined as follows: 



Xi m = -e v/2 r- l Ap lm 
Ei m = -e^ 2 r- l A Plm 
S /m - -e v l\~ l As, 



lm ■ 



(2.17) 



At this point, we already can see that the presence of the new dissipative terms affects the resulting equations in 
two main ways. First, we have a number of additional terms in the perturbed energy-stress tensor, related all of them 
to the heat flux. Second, the thermodynamical relation between the pressure and the energy density is now more 
complex because the equation of state (EOS) is not only a function of the energy density, but also of the entropy, say, 
p = p(p, s). For clarity, we will discuss the above points separately. 
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1. Taking into account the heat flux in the stress-energy tensor. 



The new terms appearing in the perturbed equations can be greatly simplified by working in a suitable reference 
frame. Introducing the four-velocity of an observer in the frame in which the energy (as opposed to the matter) is at 
rest 

= 5u » + p( o) + p (o) ' ( 2 - 18 ) 
then, the stress-energy tensor (2.10) has formally the expression of the perfect fluid case: 

ST^ = (Sp + 5p)u^u^ + 6pg$ + p^Sg^ + (Ju M u<°> + 5u v u^ ) . (2.19) 
Accordingly, we can redefine the original variables of Lindblom & Detwcilcr in the following way: 



W lr , 



W lm + e~^ , n ?V m ™ ( 2 - 2 °) 



Vim = Vl m - : l^ 2 ! m (n\ ( 2 ' 2 1) 



it = V „v/2-\ K P ( >' Qllm 

iuir pW' + p^> 



(2.22) 



and now the harmonic expansions of Su^ and of the thermodynamical quantities are also formally identical to those 
in the perfect fluid case: 

M M - f --e-^ 2 ^o/ m ^ m , iw W lrn Y lm , -\uj—^V lml ah Y l ^\ r l e^ 

(2.24) 

X lm = -e v ' 2 (r- l S Plm + e —W lm p( a A (2.25) 



Ei m = -e"/ 2 ( r- l Spim + e —W lm pW ) . (2.26) 



r 

■A/2 



Thus, we can follow the formal derivation of the perturbation equations in the adiabatic case, but using variables in 
the energy frame, the hatted variables; in this way we take easily into account the extra terms in the stress-energy 
tensor. The final system of equations (see section IV) will be analogous to the original system in the formalism of 
LD, if the perturbation of the energy E is left explicitly in the equations. 

2. Relation between Ap and Ap 

In the adiabatic case, the system of equations is closed by using the EOS to express the energy perturbation in 
terms of the pressure perturbation, i.e. Ap = c~ 2 Ap where c 2 is the adiabatic speed of sound 

This relation is used in the derivation of the LD equations and by many other authors in order to eliminate the Ap 
from the equations, so that in the final system only the pressure perturbation appears. 

For general EOSs (of the form, i.e., p = p(p, s)), the Lagrangian perturbation of the energy density can be expressed 
in terms of the Lagrangian perturbations of the other variables by using the thermodynamical relation: 

dp = c~ 2 dp + nTaids , (2.28) 

where T is the matter temperature, n is the particle density (baryon density for our purposes in neutron stars) and 
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Therefore, the Lagrangian perturbation of the energy density is given by 

Ap = C - 2 Ap + n {0) T^a 1 As. (2.29) 

In the following we will omit the super-indexes that label background quantities (n, T, p, p) for clarity. More details 
about the equation of state and thermodynamics are given in the next section. 

The perfect-fluid limit is recovered when the element of fluid does not exchange heat with its surroundings, i.e., 
As = 0. In the general case, As does not vanish, so we have to take into account all terms in equation (2.29), which 
couples the entropy perturbation to the other perturbations. Therefore, instead of simply substituting Ei m — c~ 2 Xi m , 
one must now use 

Eim = c~ 2 Xi rn + nT a x % m . (2.30) 
But this introduces a new variable, the perturbation of the entropy, so that we need additional information. 



3. Using the energy and baryon conservation equations 



The new variable which is now included in our system of equations can be rearranged by using the first law of 
thermodynamics and the continuity equation, which was not necessary in the perfect fluid case. The time component 
of the divergence of the stress-energy tensor, i.e. the energy density conservation equation, up to first order in the 
perturbations, reads 

= -^e-^ 2 Ap - (p + p)u% - (5q% + Sq a uf } u a . p ^ = . (2.31) 
In order to simplify this expression we can use the equation of continuity (2.2) 

(nu* 1 ).^ = u^n^ + nu^ = er vl2 \ujAn + nu^ = (2.32) 
and the first principle of thermodynamics in the form 

Ap = ^-^An + nTAs , (2.33) 
n 



to obtain 



that 2 , in terms of E; m , becomes 



iujnTAs = -e v l 2 (s^.j + 5q a u^u a .^ (2.34) 



iunTHim = 
,1(1 + 1) 



.1 + 1 v 1 -A' k' nTs'\ . 



Q 



2lm 



(2.35) 



Equation (2.35) gives the entropy perturbation in terms of the perturbed energy flux Q\. In order to close the 
system, we need to supplement our system with the equations that describe how the heat flux depends on the local 
gradients of the perturbations: the transport equations, which are given by the relativistic theory of dissipative 
thermodynamics. 



2 Notice that when 8q a = 0, equation (2.34) gives As = 0: the adiabaticity of the perturbations directly follows from the form 
of the perfect fluid stress-energy tensor. 
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III. A SIMPLE TWO-PARAMETER EQUATION OF STATE. 



In this section we describe the EOS that we have used and we define all the thermodynamical quantities that 
appear in the final form of the equations. We begin by considering a general equation of state depending on the 
baryon density n and the entropy per baryon s: 

p = p(n,s) ; p = p(n,s) . (3.1) 

Let us define the specific internal energy e (energy per particle, excluding the rest mass) and the speed of sound c s : 



e = rag 

n 



dp 



(3.2) 
(3.3) 



(m_e baryon mass). Next, let us assume that the EOS can be written as follows 

p(n, e) = ()-l)ne (3.4) 

with 7 constant. Using thermodynamical identities, one can show that (§f ) = -77 , so that (3.4) is equivalent to 

p(n,s) = K(s)(m B n)~* (3.5) 

where K(s) is an arbitrary function of integration that only depends on the specific entropy. If one chooses K (s) = 
Kq exp (7 — l)s, this leads to the perfect gas law p = nT (in units with the Boltzmann constant ks = 1). 

The form (3.5) allows to include finite temperature effects with a simple equation of state, just giving an explicit 
dependence of K(s) with the entropy. From the first law of thermodynamics 



dp ■■ 



the temperature can be found according to 



t - l (t 

n \ os 



p + p 



1 



dn + nTds , 



dp 



p d\nK 



The rest of thermodynamical quantities we need, in terms of K, are the following: 

P 



+ 



7-1 



a x = 



_L (L 

jK \K 

j_(dp 

nT \ ds 



-1 1 
H 



Cr> 



or 



7-1. 
7-1 



ds \ ds 



(3.6) 
(3.7) 

(3.8) 
(3.9) 
(3.10) 

(3.11) 



where C p is the specific heat at constant pressure. 

In order to study the effects of heat transport in the non-adiabatic case, we want to give a simple form of EOS which 
depends non trivially on the entropy profile and such that, consistently with the assumption that the background flux 
is negligible, it is isothermal in the relativistic sense (Te ly/ ' 2 =const.) 3 . This can be done as follows. By using the first 
law of thermodynamics (3.6) it is easy to prove that 



We stress that, as we said in the previous section, the equations we derive do apply also for background profiles not perfectly 
isothermal, as long as the background evolution timescale is longer than the dynamical timescale. 
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dT ds . . dp 

T Up p + p 

If we look for relativistic isothermal profiles, this implies that 

dp 



dT _ dv 



p + p 



where the last equality comes from the TOV equation (hydrostatic equilibrium). Therefore, we have 



ds 



ax- 



dp 



7-1 



1 



, p p + p 7 
and, substituting the explicit expression of the speed of sound, 



dp 



ds 



Defining x = lnp 1 1 ^ 7 , we finally obtain 

dx 
ds 



7-1 



1 



7-1/ P 



dp 



i + ^(f) 1/,p 



1 + 



7-1 



e x f(s) 



-In 

ds 



(3.12) 



(3.13) 



(3.14) 



(3.15) 



(3.16) 



where f(s) = K(s) 1 ^ . In principle this differential equation is not separable, but by numerical integration we can find 
the entropy profile s(p) that leads to isothermality in the relativistic sense, for any given function f(s). In particular, 
we have chosen 



/(*) = C + f3e s / a 



(3.17) 



that gives a constant specific heat C v — a. Here, C and (3 are arbitrary constants fixed to reproduce similar masses, 
radii, and temperatures as for newly born neutron stars. In particular, the temperature of the model that we analyze 
in the results section varies between (5-8) xl0 11 K 7 and C p = 0.5. 



IV. DISSIPATIVE RELATIVISTIC FLUIDS AND HEAT TRANSFER. 



The theory of standard irreversible thermodynamics was first extended from Newtonian to relativistic by Eckart 
in 1940. This theory shares with its Newtonian counterpart the problem that perturbations propagate at infinite 
speeds, which results in unstable equilibrium states. We must notice, however, that the non-causal Newtonian 
theory has been used for years in the study of non-relativistic stars with remarkable success, and with observational 
confirmation. Therefore, there is no reason to believe that the standard theory cannot be used as a first order approach 
to the problem and can give reasonable estimates of the main effects, even for relativistic stars. In general relativity, 
causality can be restored within the framework of an extended theory developed by Israel and Stewart [17]. This 
extended theory is also known as causal thermodynamics or second order thermodynamics because of the appearance 
of second order terms of the dissipative variables in the entropy. The problems associated to non-causal heat transfer 
become more relevant when the mean free path becomes larger than the spatial scale of the problem (in our case, 
the radius of the star), that is, for large thermal conductivities or diffusivity. The extended theory automatically 
incorporates transient phenomena on the timescale of the mean free path that cure the inconsistency. In this paper, 
we are interested in understanding qualitatively the effects of dissipative terms on the oscillation properties of stars, 
rather than developing a fully consistent theory, which is a more ambitious and complex task. Nevertheless, we start 
with the extended theory of dissipative fluids, and later on we discuss which second order terms are neglected and 
their relative importance. 

We begin with the equation that describes the variation of the entropy 

nTu a s, a = -q a . a -q a uf 3 u a . :0 , (4.1) 

and the definition of the heat flux, including the relaxation term that restores causality; this is of covariant Maxwell- 
Cattaneo form (a truncated version of Israel-Stewart equations): 
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nh^q^ + q a = -K{h%r if) + TvPu a .£ 



(4.2) 



where h@ is the projector onto the spatial subspace, h@ = 5@ + Uavfi , k is the thermal conductivity, and T\ is the 
relaxation timescale in which the system restores thermal equilibrium. Let us consider equations (4.1), (4.2) on a 
perturbed static spherical background. The perturbation of equation (4.1) gives the equation previously derived 
(2.35). The perturbed version of equation (4.2) must satisfy q a u a = and, since we are imposing a zero background 



flux (q^ =0, j = 1, 2, 3), we have 5qo = and 

(1 + iwe-^ 2 Ti)(% = — k (STj + STu^ufl + T^S(u^ Uj ,^ . 
Expanding in spherical harmonics and after some manipulations we find 



(1 + iuer l//2 T 1 )5q J = -e 



5 (Te v ' 2 ) 
A (Te^ 2 ) 



- KTcj 2 e~ v Vj 



KTuu 2 e- v v j 



(4.3) 



(4.4) 



where 



e A/2 

rH Um Wi m 



Y lm ,V lm YT )r l e 



Applying the thermodynamical relation (3.12), which implies 

AT As A P 
-7jr = TT + (l-«i)- 



P + P 



(4.5) 



(4.6) 



on Eq. (4.4), changing to the energy frame (perturbations in terms of hatted variables), and applying the Einstein 
equation 



X lm -u 2 (p + p)e-^ 2 V lm - ^-e^/ 2 (p + p)W lm - l(p + p)e^ 2 H 0lm = 0, 



(4.7) 



we obtain a first algebraic relation between the entropy perturbation, Q2, and X, 



(l + iu J e-^ 2 T 1 )Q 2l 



-Xim ■ 



C p p + p p + p 



12 lm 



(4.8) 



The exact evaluation of the relaxation time T\ involves complicated collision integrals and depends on the microscopical 
details of the interaction between particles. Notice though, that it is usually estimated as ~ kT ' j p (see i.e. [17]). 
Therefore we do know that this second order correction is of the same order as the last term on the right hand side 
ofEq. (4.8). Indeed, if one defines the relaxation timescale as n —nT/{p + p), these two terms cancel each other. In 
the following, we will neglect all second order terms, and Equation 4.8 becomes 



Q 



2lm 



Cn 



ai 



X, 



I in 



(4.9) 



We have checked that the numerical results presented in section VI are not altered by including the second order 
corrections, being the relative difference between the two solutions less than 1%. 

Let us come back to Eq. (4.4). After neglecting the second order term, the quantity e u ^ 2 Sqj/ ' k + Tbj 2 e~ u t 2 Vj is a 
gradient, thus the harmonic components of the perturbed heat flux satisfy 



Q2 



2lm 



'-Q2lm + -Q21 

r r 



-Tco 2 e-^ 2 



-)V lm - rH 

r 



e V2 „ 
Urn H W h 



(4.10) 



Next, we can make use of the following equation, that is derived directly from the LD equations [12,13] 
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Wlm Vh, 
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-v/2 



{p + p)' 



:{Xl m p' - Elmp') 



to simplify the right hand side of Eq. (4.10) 
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Q'2lm = Qllm + -Qllm + 7 j 7^(Xl m p' - Ei m p') 

r r \P + P) 



(4.11) 



(4.12) 



and, finally, introducing the thermodynamical derivatives defined in the previous section (ot\,C p ) and keeping only 
first order terms we obtain 



Q'llra = --Qllm. + -Qllm + nTa\ — - ^rQ2lm ■ 

r r 2 (p + p) 



(4.13) 



V. NUMERICAL IMPLEMENTATION 



A. The complete set of equations 



We can now write the full set of equations. For simplicity, in the following we will omit the sub-indexes Im, keeping 
in mind that one has a complete set of equations for each multipolc in the expansion. Our variables are 



H , H u K, W, V, X, E, S, Q u Q 2 . 
Four of these variables Hq, V, E, S, are given in terms of the others by the following algebraic relations: 



1 



3M + C-i)('+2) r + 47rr 3 p 



»Kr 3 e- v/2 X - r 3 e- x 
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(i-m+2) 



,2^3„-i/ 



rV 



r — u>~r"e ' — e " | — 1 

2 2 



rv 



2 2r 
K 



Hi 



(5-1) 



(5.2) 



X e -u/2 e -A/2^ 1 
JpTp) ~—J W -2 Ho 



E = c- 2 X + nTaiS 
'ft "l 



S — C„ 



T p + p 



X 



(5.3) 
(5.4) 

(5.5) 



The remaining variables H\, K, W, X, Qi, Q 2 , satisfy the following differential equations: 



1 

r 

r 



1 + 1 + 2M— + 4irr 2 e x (p - p) 



Hi 



Hq + K - 16n(p + p)V 



r 2r 



i + i i ; 
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(5.6) 
(5.7) 
(5.8) 
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Q 1+ e xl l±ll 



(5.9) 
(5.10) 
(5.11) 



It must be noticed that the first four equations are formally identical to those of Lindblom & Detweiler [12,13], but 
in terms of the new variables defined by Eq. (2.20)-(2.23). 



B. Boundary conditions 



We seek for solutions of the perturbation equations which are regular at the origin. Assuming that all variables 
near the center of the star have the form x{r) = x(0) + 0(r 2 ), and expanding equations (5.2)— (5.11) we find that the 
following relations must be satisfied at the origin: 



H (0) 
V(0) 



K(0) 
-jW(0) 



X(0) = ( Po +Po)e 



v /2 



4tt, „ . 2 e- U0 
y(po + 3p ) -w — — 



W(0) + -K(0) 



Hi(0) 

Q 2 (o) 



1 



1(1 + 1) 

yQi(O). 



21K(0) + 16tt( Po +Po)W(0) 



(5.12) 
(5.13) 

(5.14) 
(5.15) 

(5.16) 



The values of Ei m and S; m at the origin are given by (5.4), (5.5) evaluated at r = 0. Therefore, our equations admit 
three independent solutions, which are reduced to one by imposing appropriate boundary conditions at the surface of 
the star. The first is the vanishing of the Lagrangian pressure perturbation 



X(R S ) =0: 



and the second is the vanishing of the radial flux 



Qi(R s 



0. 



(5.17) 



(5.18) 



Once we know the only independent solution at the surface of the star, the Zerilli function and its derivative can be 
computed in terms of the Hi and K, as usual (see for example [18] and [19] for the expression of the Zerilli function 
in terms of H 0} K; the expression in terms of Hi follows by applying Eq.(5.2) ): 



Z = r l 
Z' = r l 



2r 2 



:(K-e»Hi) 



+e 



(I - l)(/ + 2)r + 6M 
2(1 -!)(/ + 2)r 2 - 6M((Z -!)(/ + 2)r + 2M) 
((l-l)(l + 2)r + 6M) 2 
2 „ (I - 1)1(1 + l)(l + 2)r 2 + 6M((l + 2)r + AM) 



K 



((I - l)(l + 2)r + QMf 



(5.19) 



Notice that the metric variables Ho, Hi, K had not been redefined by the presence of heat flux. Finally, we integrate 
the Zerilli equation outside the star using the continued fraction method (see e.g. [20,21]), to obtain the amplitude of 
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the ingoing and outgoing parts of the gravitational wave. The quasi-normal modes will correspond to the solutions 
for which the wave is purely outgoing. 
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FIG. 1. Representation, in the complex plane, of the variation of the frequency and inverse damping time (imaginary part 
of the frequency) as the coefficient n/n varies from the adiabatic limit (<C R) to the isothermal limit (> R). The results 
correspond to the first p— mode. As expected, the thermal damping is maximum when n/n is of the order of unity. 



VI. RESULTS 



In the previous sections we have ignored the effects of chemical diffusion. Our equations can then be applied to 
a physical situation in which heat conduction is relevant, while chemical diffusion is not. This could be the case 
of a young, hot, neutron star about 30 seconds after birth, once the proto-neutron star has lost its lepton content, 
and the temperature gradient is not far from isothermal. It is also the situation that one expects to find in a cold 
neutron star that undergoes a phase transition to deconfined quark matter in a color superconductor state. In this 
latter scenario, pairing between quarks inhibits the presence of leptons [22]. For reference, in Appendix B we derive 
the equations including the terms related to chemical diffusion, leaving for future work their application in a realistic 
scenario. In this paper, we will focus only on the effects of finite thermal conductivity, from which we can understand 
most of the qualitative differences between the adiabatic and non-adiabatic cases. Hereafter we consider a profile for 
the thermal conductivity such that n/n is constant throughout the star. The ratio n/n represents, roughly, a sort of 
mean free path, to be compared with the characteristic scale of the system, i.e., the radius of the star. In the diffusion 
limit n/n <C R, perturbations can be considered basically adiabatic, while when n/n becomes of the order or larger 
than the radius of the star, thermal equilibrium is achieved in a timescale shorter than the dynamical characteristic 
timescale. 

We begin by studying the effects of finite thermal conductivity on the p— modes. In Fig. 1, we show, in the complex 
plane, the variation of the complex frequency as the ratio n/n varies from the adiabatic limit to the isothermal limit. 
The oscillation frequency (real part) is represented in the horizontal axis while the vertical axis corresponds to the 
inverse damping time (imaginary part). From the results, we can see that the real part of the frequency of the first 
acoustic mode is shifted to lower values in about 200 Hz, as the mean free path becomes larger, an the isothermal 
limit is approached. The damping time in both, the adiabatic and isothermal limit is very similar (not visible in the 
scale of the figure), being in both cases of about 1 s. However, in the semi-transparent regime, thermal dissipation is 
so effective that the damping time is reduced by about 3 orders of magnitude (to 2-3 ms). These results can be easily 
understood by looking at a simple toy model that describes how acoustic modes are affected by thermal diffusion. 
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We describe this simple Newtonian problem in Appendix A. It makes more manifest the relevant physics without the 
complications of the full set of equations in General Relativity. By analogy with this toy model, one can understand 
that the oscillation frequency in the isothermal limit is smaller than that of the adiabatic case, depending on the ratio 
between the adiabatic speed of sound c s and the isothermal speed of sound Ct, defined as 



1 nTaxCpioix - 1) 

C 2 s P + P 



(6.1) 



For the EOS we used, this difference is about a 5% in average, which is consistent with the shift in the p— mode 
frequency. 

Several comments about the adiabatic and the isothermal limits are in order. In both cases, the heat flux tends 
to zero as we approach the limits, but for different reasons, that can be understood by looking at Eq. (4.9). In the 
diffusion limit (adiabatic perturbations) k — > 0, but in the isothermal limit, the flux vanishes because the term within 
brackets (i.e. the perturbation of the temperature, including relativistic corrections) vanishes. Therefore, in both 
limits Eqs. (5.10)— (5.11) are decoupled from the rest of the system. The difference between both limits appears only 
in Eq. (5.8). Explicitly, in the adiabatic limit Eq. (5.8) reads: 



W 1 = -—W + re x ' 2 



p + p cj r z 2 



(6.2) 



while in the isothermal limit it becomes 



w' = -i±lw 



re 



A/2 



er v ' 2 X _ nTaiCp Av _ IJJ_ 
p + p c\ p + p 2 i 



(6.3) 



where 



e -vv - 

Av = H Q + W. 



(6.4) 



In other words, one just needs to substitute the adiabatic speed of sound by the isothermal speed of sound (consistently 
with what happens in Newtonian perturbation theory) but now there appears an additional term (proportional to the 
Lagrangian perturbation of v/2) because of the fact that the concept of isothermality in general relativity involves 
the red-shifted temperature, Te u ^ 2 . 

The damping time when dissipation is most effective, i.e. (as we can see from Fig.l) when n/n w 1 km, can be also 
estimated by 



Tdiss — 



ds/dt 



= nC, 



R 



(6.5) 



and taking typical values at the interior of the star, such as C p = 1, R = 10 km, we get that for n/n — 1 km, the 
dissipative timescale is Tdiss = 0.001 s, in agreement with the results shown in Fig. 1. These results show explicitly 
and quantitatively that thermal dissipation affects the damping of the non-radial pulsations, competing with the 
other main dissipation mechanism, i.e., GW emission. The damping times we compute take for the first time into 
account both mechanisms in a self-consistent manner, by including the appropriate physics in the equations. In the 
semitransparent regime, when diffusion is more effective, we must expect short-lived, strongly damped, GW signals, 
more similar to a GW burst than a proper ~ kHz oscillation lasting for several hundred oscillations. Notice, however, 
that if the physical conditions are such that the thermal relaxation timescale becomes shorter than the oscillation 
period, thermal damping does not affect much the damping time, being the only remarkable difference a shift in the 
oscillation frequency. 
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FIG. 2. Same as Fig. 1 but for the fundamental mode. 

Let us now focus on the fundamental mode. The /—mode is a global oscillation mode that depends essentially 
on the average density of the star, rather than on the velocity at which acoustic waves propagate. In fact, it is well 
known that an incompressible fluid (infinite speed of sound) does not have p— modes but the frequency of the /—mode 
is similar to that of a realistic star with the same average density. For these reasons, one expects similar results to 
those of the p— mode but with smaller effects on the shifts of frequencies. In Fig. 2 we show our results for the 
/—mode, when this is clearly visible. Now, the shift in frequency is only of a few Hz, and the effect on the damping 
time is also smaller, although in the semitransparent regime, it can be as much as a factor 2.5 smaller than that of 
the adiabatic case (0.25 s). Analogously to what happened for the p— mode, once we enter in the isothermal regime, 
thermal damping is less effective and the damping time becomes quite similar to the adiabatic value. One must keep 
in mind that these results depend on the particular details of the EOS, or more precisely, on the temperature and 
specific heat. In cases with higher temperatures and lower specific heat, the frequency of the fundamental mode could 
be quite different from the adiabatic case. It remains to be analyzed what are the quantitative differences for realistic 
models at different stages of neutron star evolution, which is beyond the scope of this paper. 
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FIG. 3. Same as Fig. 1 but for the first g— mode. 

Finally, let us discuss the g— modes. In Fig. 3 we show the results for the first gravity mode, which in the adiabatic 
case has a frequency of 364 Hz and a damping time of thousands of seconds. The existence of these modes is related 
to the presence of thermal (or chemical) gradients in the star, therefore it is not surprising that as we increase 
the thermal conductivity and, consequently, heat interchange between displaced fluid elements is more effective, the 
g— mode frequencies decrease and the damping time increases. Contrary to what happens with the acoustic and 
fundamental modes, in the limit of infinite thermal conductivity the g— modes degenerate at zero frequency. It is 
interesting to note that the damping time does not decrease indefinitely as the mean free path increases, being 
bounded by the minimum dissipativc timescale that is, again, close to the estimate of Eq. 6.5 (of the order of ms), in 
our scenario. The last important difference with the other modes is that for values of n/n as low as 10~ 3 km, although 
the change in frequency is small, the damping time is sensibly shorter (0.01 s). For the /— or p— modes it was needed 
to have higher values of n/n to obtain significant differences with respect to the adiabatic damping time. 



VII. FINAL REMARKS 



In this paper we have developed the formalism to study non-radial oscillations of rclativistic stars in the frequency 
domain giving one step forward from the perfect fluid case by including the effects of thermal diffusion in a fully 
relativistic formalism. This allows us to understand one of the non-adiabatic processes (we did not consider viscosity) 
that may be relevant in the study of the oscillation properties of newly born neutron stars, strange stars, or neutron 
stars with deconfined quark matter in the core (hybrid stars). When a general equation of state that depends on 
temperature is used, the perturbations of the fluid result in perturbations of temperature (or chemical composition) 
and, consequently, in heat flux (or chemical diffusion), that is coupled with the geometry through the Einstein field 
equations. We have analyzed separately the p— modes, the fundamental mode and the g— modes, each one being 
affected in a different way by thermal diffusion. As expected from a simplified model discussed in Appendix A, the 
p— mode frequency is shifted to lower values in a factor roughly proportional to the ratio between the isothermal 
and adiabatic speeds of sound ct/c s1 thus depending very much on the EOS employed and the values of the local 
temperature. The frequency of the /—mode, however, is barely affected because it is a global oscillation mode 
that does not depend much on the local value of the speed of sound. Both /— and p— modes are more efficiently 
damped in the semitransparcnt regime, when the mean free path of the particles responsible of the heat transfer is 
smaller, but close to the typical length scale of the system. The reason is that once the timescale to reach thermal 
equilibrium becomes shorter than the oscillation period, the fluid oscillates keeping the temperature constant, and 
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further increasing the thermal conductivity does not change this situation. In this limit, there is no additional thermal 
damping at first order (as shown in Appendix A), so that the damping times in the adiabatic and isothermal limits are 
very similar. The response of the g— modes to heat transfer is quite different. Since these modes exist because of the 
presence of thermal gradients, it is naturally found that, as heat transfer becomes more effective, and the temperature 
perturbations are smeared out, the frequency is shifted to lower values and the damping time becomes shorter. In 
the isothermal limit the g— modes have degenerated to zero frequency. 

In this first approach to the problem, we have only presented results about the effects of heat transfer. It must 
be remarked that in newly born neutron stars the effects of lepton diffusion are also important, and in Appendix B 
we discuss the relevant equations to study that case. The coupling between thermal and chemical diffusion is one of 
the complications that arise if one wants to study the realistic case of proto-neutron stars a few seconds after birth, 
but there is a second one. In this early stage the background gradients of temperature and chemical potential are 
important, and in some cases, or for some of the modes, the overall evolution timescale may become similar to the 
characteristic oscillation periods or damping times. Furthermore, proto-neutron stars are expected to rotate fastly, 
while in the present work we focused on thermal effects neglecting rotation. We defer for future work a more rigorous 
study of the realistic scenario, that must consider the non-trivial issue of coupling between rotational and thermal 
effects. This must probably be done in the time domain instead of the frequency domain, or including second order 
terms. 

Our results can also be directly applied to another interesting case: hot strange stars in which quark matter is in 
a color superconductor state. As discussed in the paper, this appears to be the natural state of deconfmed quark 
matter, and recent calculations show that the thermal conductivities are many orders of magnitude larger than those 
of standard neutron star matter. In this latter scenario, if the strange star is born after the mini-collapse of an old, 
evolved neutron star that has been accreting matter, the initial lepton content is small, and the effects of thermal 
diffusion are the dominant non-adiabatic correction. 
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APPENDIX A: A TOY MODEL FOR NON-ADIABATIC OSCILLATIONS. 



Consider a 1-dimensional problem consisting of a box at constant density, and constant pressure, and let us study 
the normal acoustic modes of the fluid. The linearized continuity and momentum equations in the adiabatic case are 

^+pVt* = , ^T + -V P1 = 0. (Al) 
at at p 

Integrating in time the continuity equation gives pi + pV£ = 0, with £ = Sr, and substituting pi in the momentum 
equation leads to 

^ = —V P i = -£vpi = & 2 Z (A2) 
dt z p p 

which is a simple wave equation that, when we consider perturbations of the form exp(iat — ikx) gives the dispersion 
equation: 

-a 2 + k 2 c 2 s = 0. (A3) 

Its solution consists of oscillatory modes with frequencies ±kc s . 

Consider now the non-adiabatic case. The perturbation of the pressure is 

pi = c 2 Pl + pp Sl (A4) 

with /3 = j(dp/ds) p , and equation (A2) becomes: 



15 



g = c 2 V 2 £ /3V S1 . (A5) 

We have to introduce the additional equation of conservation of energy 

ds 

nT— = -VF (A6) 
at 

where T is the temperature, n = p/m is the particle density, and 

F = -kVT. (A7) 
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FIG. 4. Left panel: Solution of the dispersion equation (A13) as a function of D = nk/nCv-The oscillation frequency, 
normalized to kc s is represented by the solid line, while the inverse damping time (imaginary part of the frequency) is denoted 
by 1/t (dashed line). Right panel: representation in the complex plane of the real (frequency) and imaginary (inverse damping) 
parts. 



Considering first order perturbations of the previous equations we have 

^dsi „~ 
at 

that can be used to replace after taking the time derivative of Eq. (A5), to obtain 



?|=^V 2 ^-^V(V 2 T l) . 
dt 3 s dt nT ' 



(A8) 



(A9) 



Next, we only need to find an expression for the last term in Eq. (A9). The perturbation of the temperature can be 
written in terms of p\ and si as follows 



Ti 



dT\ 2 T 
dp J . c v 



where c v = T (■§§;) p is the specific heat at constant volume. Thus , 



dT 



c 2 Vpi + — Vsi = - 



Inserting this in Eq. (A9) one obtains 



/dT 
\~d~P~ 



*1 

I3c v \dt 2 



d 3 ^ 
dt 3 



dt nT \ dp ) , 



K 

nc v 



ft 
dt 2 



-c 2 V 2 C 



(A10) 



(All) 



(A12) 



1G 



Considering again perturbations of the form exp(iat — ikx) we can derive the dispersion equation: 



-a 3 + c 2 k 2 a + i— k 2 \a 2 -c 2 T k 2 } =0. (A13) 
8 nc v 1 1 1 

The adiabatic limit is recovered when n/nc v -C 1/k while in the limit k/uc v 3> 1/fc, the dispersion equation becomes 

cr 2 - k 2 c 2 T = , (A14) 

which simply states that sound waves propagate now at Ct instead of c s . It must be remarked that in the isothermal 
limit, the solution is always real (no damping), and a acquires an imaginary part only in the intermediate regime. 
Notice also that dissipation is most effective for modes at short wavelengths (large k). In Fig. 4 we show the solutions 
of the dispersion equation (A13) as a function of the parameter D — nk/nc v . The oscillation frequency (normalized 
to kc s ) is represented by the solid line, and the inverse damping time by the dashed line. The qualitative behaviour is 
very similar to the results discussed in the text. The ratio of the the limiting frequencies depends on the ratio ct/c s . 
For this example, we have taken ct/c s — 0.8. 

A final remark is in order. The dispersion relation can be written as 

7 3 + «27 2 + ai7i + «o = (A15) 

a 2 = D , ai = c 2 , a = Dc 2 ^ (A16) 

where 7 = ia/k. The condition that at least one of the roots has a positive real part (unstable) is equivalent to one 
of the following inequalities 

a 2 < 0, a < 0, aia 2 < a . (A17) 
Since all coefficients are positive defined, only the third condition could apply. It is simply 

c 2 < 4 (A18) 
which is generally not the case, for reasonable realistic EOSs. 

APPENDIX B: INTRODUCING NEUTRINO DIFFUSION IN OUR EQUATIONS 

In order to study non-adiabatic perturbations in proto-neutron stars, we need to generalize our equations to the 
case when both, energy flux and particle number flux due to neutrino diffusion are present. In this appendix we derive 
the equations in a similar way as it was done in the main body of the paper for the case with only heat transfer. 

1. Thermodynamical variables and relations 

We choose as independent thermodynamical variables the pressure p, the entropy per baryon s, and the lepton 
fraction Yj, — n-^jn, with n being baryon number density and the lepton number density. Given these three 
thermodynamical quantities, the EOS gives the rest of variables. In particular, the EOS provides 

n = n(p,s,Y L ) p = p(p,s,Y L ). (Bl) 

The reason for taking Yl as the extensive variable describing the chemical composition is the following. The work 
associated with a variation of chemical composition at constant entropy is 

Y^IMdYi. (B2) 

i 

Here, Yi — rii/n, and \ii,ni are the chemical potentials and number densities of the specie i, respectively. In our case, 
the process that changes the chemical composition is the capture of electrons by protons and its reciprocal (inverse /? 
decay) 

e + p <-> n + v , (B3) 
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so that chemical equilibrium results in p, v = p, e + p, p — p n . Since only neutrinos can diffuse throughout the star, we 
also have dY e = dY p = —dY ni while dY v is independent from the others. Consequently, 

^ [lidYi = (p e + Hp - Hn)dY e + [LydYy = p, v {dY e + dY v ) = p u dY L . (B4) 

i 

Notice that in this system the variable conjugate to Yl is p v . 

According to this result, the first principle of thermodynamics (see for example [14]) takes the form 

dp = ^——^-dn + nTds + np v dYL ■ (B5) 

Analogously, since we use (p, s, Yl) as independent variables, we can write, in place of (2.28), 

dp = Cg 2 dp + nTaids + np u a 2 dYL (B6) 
where the sound speed c s is defined by 

■ (B7) 



dp 



and we have defined the following thermodynamical derivatives: 



-^(IL- -sH&i ■ <B8) 



We also define 



P,Y L 



where the last equality can be easily proved from Maxwell relations. 
Let us now consider the heat function SQ 



SQ = Tds + Y^ VidYi = Tds + p u dY L . (BIO) 

i 

Then, we define 

(§) =K^L,^ P 



p ,y l \ WJ - y p ,y l 



P ' 

p,s 



where C p is the specific heat at constant pressure and constant number fraction, A p is a coefficient describing the heat 
associated to a change in composition at constant entropy and pressure. With all the above definitions, and using the 
first law of thermodynamics, the expressions of dT and dp v in terms of s,p, Yl can be written as 

dT= — ds+ (1 -a 1 )dp + a i dY L (B12) 

C p P + P 

dp u = ^dY L + -^(l-a 2 )dp + a 3 ds. (B13) 
A p p + p 

Equations (B6), (B12), and (B13) allow us to express the perturbations of density, temperature and chemical potentials 
in terms of the independent variables. Similar relations apply for the gradients of the background quantities. For 
example, from (B6) 

Ap = c^ 2 Ap + n(°¥°)aiA« + n^p^a 2 AY L (B14) 

p (0), = c -2 p (0)/ + nT (0) aiS (0)> + nfl (0) a2 y(0), (B15) 

and the same applies for AT, Afj, u , T^', p^' . 
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2. Dissipative thermodynamics with energy and number transport 

Next, we need to generalize Eckart dissipative thermodynamics to the case when also chemical diffusion, (lepton 
diffusion in this case) is present. Therefore, we need to consider the lepton conservation equation 

{nY e u a +nY v (u a + n a )). a = 0, (B16) 

where n a is the neutrino four-velocity with respect to the frame comoving with matter (see [23], [17]), with n^u^ = 
and n M n M = 1. Defining the neutrino flux as f a = nY v n a , that satisfies u a f a = , and using the continuity equation 
(2.2) we obtain the equation of lepton conservation nu a (YL), a = —f^.p- This latter equation is particularized for 
lepton diffusion, which is relevant for nascent neutron stars, but it is a general conservation equation for massless 
particles. Therefore, the set of equations we have is the following: 

• Continuity equation: 

(nu a ). a ^0. (B17) 

• Energy conservation equation (stress-energy tensor projected onto u): 

u a p, a + u a . a (p + p)+ q a . a + q a uf 3 u a , fj = . (B18) 



• First law of thermodynamics: 



dp = ^-^-dn + nTds + np u dY L . (B19) 



Lepton conservation: 

nu a Y L , a = -f% . (B20) 



Together they give: 

nTu a s, a = -q a [a + n v f« a - q a u f3 u a[ p . (B21) 
Let us now define the entropy flux as in [24] , 

S a = snu a + ^-p^. (B22) 

Notice that when degenerate neutrinos dominate the transport of energy and particles, q a = fi v f a and there is no 
entropy flux. By taking the divergence of the entropy flux, we have 

TS% = -^(D a T + Tv?u a .fi) - f a TD a n , (B23) 

where we have defined r\ = The second law of thermodynamics, S a . a > 0, must be always satisfied. The simplest 
assumption for the form of the energy and lepton fluxes that satisfies S a a > is 

la = -ke(D u T + Tu l3 u a; p) , f a = -K N D a r) (B24) 

with ke, kn positive defined coefficients, which are, respectively, the thermal conductivity and the diffusivity governing 
the transport of particles. In the static, radially symmetric case, our equations coincide with a particular case of those 
of [8]. 
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3. Perturbations on a static, radially symmetric, non stratified background 



Let us consider a static, non-stratified background space-time, in which the lepton and energy fluxes satisfy 

/i 0) =0 g£ 0) = . (B25) 

Consistently, we will have (Te"/ 2 )' = 0, rj' = 0, (p v e v / 2 )' = 0. Next, we expand the fluxes in spherical harmonics as 
follows: 



5q. t = K E e-^(Q Um r l - L Y lm ,Q 2lm r l Y^)e iut 

„2 



Sfi = un- 



it 



(F Um r l - 1 Y lm ,F 2lm rXTy wt , 



l\rlm\Aujt 



as well as the Lagrangian perturbation of Yl 

AY L = -e- v ' 2 r l y lm Y lm e 10jt . 
Analogously to the procedure in subsection 2.1, we define 

-0 -., „u/2-\ k eY[ Qum 
yim = yim - e — ; 

iuir p + p 

By expanding the transport equations (B24), we obtain 

"A(Te"/ 2 ) 



(B26) 



(B27) 



(B28) 



% = -e-"l 2 TK E 
S ft = -K N Arn , 



K E Tu> z e~ v Vi 



(B29) 



where Vj has been defined in (4.5). 

At this point, we can derive the equations for the perturbations E, Qi, Q 2 , F\, F 2 , S, y. The differences with respect 
to the case with only thermal diffusion are the following: 



• Expression for E. 
From (B6) we have 



Eim = c s X im + nTa\T,i m + n\i v a 2 yi r . 



(B30) 



• Equations for Qi, F\. 
Perturbing the conservation equations (B21), (B20) we have 



nTAs + np, v AY L 
nAY L 



e v ' 2 

1U) 



(dq a a + 5q a uPu a , ) 



that, after expanding in harmonics, give 

nT'Eim + np v yi m = k e — 

\ur 

x l(l + l) 



—e 



Ql lm 
Q2lm 



l + l v' — A' k' e nTs' + np v Y' L ' 

Qllm 



K E 



P+P 



np v yi m = K N r) — 



ie v/2-X 



,1(1 + 1) 



F[lr, 



l + l V>-\' K' N 
1 ~ 1 #llm 

r 2 k n . 



np v Y' L Qi k e 
p + p rj 2 k n 



(B31) 
(B32) 



(B33) 



(B34) 
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• Expressions for E, y. 
Comparing the i = 9,(p components of (B26) and (B29) we obtain 



-ullQllm . 2 -i/t> ( ^ T lm . Ai// m \ 

e ' — ^-+u z e "Vj^ = - I — — + J 



-"/2_1 
Pu 



(B35) 



Then, using the thermodynamical relations derived in Section B 1, the Einstein equation (4.7), and the fact that 
the background is static, that leads to 



— - -^_n' + —Y' - 

c p p + p p + t Yl -° 



(B36) 



we find 



E( m H — %^r-yim = C P A E 

Sim H — -E; TO = A„An 

Pv 



(B37) 
(B38) 



where 



P + P 

The solution of this system is 



^i-E — ; Mm H T^7~ 



An — — ^—Xir. 

P + P 



E; m — 



yirn = 



1 ~ Tfj,„ 



A„ 



7>„ 



A E - 



^2 lm + Q2 lm 



i-N 



osC, Ae 

Pv 



(B39) 



(B40) 



• Equations for Q2,F 2 . 

The equations for Q 2 , F 2 follow from the fact that e v / 2 5qj/n E + Tuj 2 e~ v / 2 Vj and 5fj/n N are gradients, which 
can be used to relate the radial and angular components of (B26), 



I 1 

^2 lm ~ F2 lm H "l Zm • 

r r 

, Z 1 T - , - , 

( ?2im — Q2im + ~Q\lm + 7 | 7^(MmP ~ El m p ) . 

r r (P + P) 



(B41) 
(B42) 



The last term in (B42) can be expanded in terms of the rest of variables. Using (B15) and (B30), and after 
some algebra, one finds 



Q f 2lm — Q2lm + ~Qlh 

r r 
v> 1 



2(P + P) l 
+np ly a 2 A p 



7>„ 

_ a 3 C p 



Pv 



nTa\C p 



Q2lm+F : 



a 3 Ar 



n 12K p tp 

W2lm T^~^2lm 



2 h. 



(B43) 
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Equations (B30), (B33), (B34), (B37), (B38), (B41), (B43), together with equations (5.2) to (5.9), form a closed 
system that allows to compute the non-radial perturbations of a star with energy and lcpton transport. 
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